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Abstract 

Autoresonance is a phase locking phenomenon occurring in nonlin- 
ear oscillatory system, which is forced by oscillating perturbation. Many 
physical applications of the autoresonance are known in nonlinear physics. 
The essence of the phenomenon is that the nonlinear oscillator selfadjusts 
to the varying external conditions so that it remains in resonance with 
the driver for a long time. This long time resonance leads to a strong 
increase in the response amplitude under weak driving perturbation. An 
analytic treatment of a simple mathematical model is done here by means 
of asymptotic analysis using a small driving parameter. The main result 
is finding threshold for entering the autoresonance. 



1 Introduction 

We consider a simple mathematical model of forcing oscillations given by the 
nonlinear ordinary differential equation 

u" + (1 + 7M 2 )u = 2af(t) cos(<p(*)), t > 0; < a < 1 (1) 

where the right hand side represents a small external force. The zero initial 
condition is added here 

(«,u')l*=o = (0.0), (2) 

so the system is starting from stable equilibrium. A nonlinearity may have a 
different signs: 7 = ±1. The driving amplitude / is a slow varying function in 
contrast to the phase function so that f'/ip'(t) — 0(1), a — ► 0. In the paper 

*The russian version of the paper is presented in "Doklady Akademii Nauk" ("Doklady 
Mathematics" in English) 
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we derive the criterion of the autoresonance phenomenon. Namely, we find a 
condition under which the system's energy £(t, a) = [(u') 2 + u 2 }/2 + 7u 4 /4 — 
u2af cos(ip) grows up to the order of 0(1) as t — > oo while the driver is being 
small: < a < 1, f (t) = 0(1). 

Great growth in the energy (and in the swing amplitude as well) usually 
takes place due to the resonance in oscillating system. The simplest model of 
the resonance phenomenon is given by the linear equation of harmonic oscillator 
under periodic force 

u" + uj 2 u = cos(i4). 

The energy of this system grows infinitely in time, if the driving frequency 
coincides with the free frequency uj 2 = v 2 . In the case of nonlinear system 
such type of phenomena is referred to as autoresonance [1-4]. Unlike linear 
systems nonlinear systems must have two distinctive features to be able to enter 
autoresonance. The first one is that the driving frequency must be varying, 
because the free frequency, depending on the energy, is varying in accordance 
with the growth of the oscillating amplitude. The second feature is that the 
driving amplitude must exceed a sharp threshold, depending on the chirp rate 
of the forcing, as it was noticed recently [1,2]. 

An analytic treatment of the threshold phenomenon is done here by means 
of asymptotic analysis of the problem (1),(2), using a small parameter a. We 
believe that an accurate study of the simple example provides understanding of 
more complicated problems. In particular, the fact of the cubic nonlinearity in 
(1) does not play any role. Similar results take place for each oscillating system 
with a smooth nonlinearity. 

As was pointed out above the free frequency depends on the energy in non- 
linear system. In order to reach large energy of the swing, using a resonance 
effect, it is necessary to vary the driving frequency d<p/dt in time, starting from 
initial value dip/dt\ t=0 = 1. In this way one has to take into account that energy 
is small and varying slowly through initial stage until the amplitude, starting 
from zero, increases up to the order of unity. Hence there are two additional 
small parameters in the setting of the problem. The first one is the character- 
istic frequency mismatch, the second is the chirp rate, i.e. the rate of change 
of the mismatch. We shall identify both small parameters as e and e x in the 
phase function: 

tp = t + £- A $(£ 1+A t), (0 < e < 1, A = const > 0), (3) 

In this work the small parameter e will be related to the amplitude parameter 
a and two types of the data /, ip(t, a) will be considered. Two different cases 
are distinguished one from another by the rate of change of the slow varying 
functions /, ip(t, a) and they are sometimes called as "the rigid frequency chirp- 
ing" and "the loose frequency chirping" . In each case a threshold for the driving 
amplitude is found and entering the autoresonance is proved in the paper. 

We consider here just only the initial stage of amplitude rising and omit the 
following stage at which the evolution of the large amplitude of the order of 
unity occurs. 
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Our theory may be considered as an asymptotic analysis of a small amplitude 
solution. This approach gives an appropriate tool in finding the autoresonance 
phenomena in cases under consideration. Note, similar results for both prob- 
lems of the nonlinear resonance and of the synchronization of oscillation [5-7] 
are known. However, they are not directly related to the autoresonance phe- 
nomenon. 



2 Anzatz of asymptotic solution 

Our approach is a version of two scale methods, [8-13]. The WKB type anzatz 
is taken as an asymptotic solution of the problem (l)-(3): 



u = a 1 ' 3 



A(et,s)exp(itp(t,e)) + c.c. + a 2/3 u 1 {t,e) + 0{a 4/3 )}, a -> 0. (4) 



This is a small amplitude approximation because of the factor a 1 / 3 . The expo- 
nent 1/3 is chosen here in order to take into account both the nonlinearity and 
the driver in the slow modulation of the amplitude A. For similar reasons in 
order to take into account the frequency mismatch <p' — 1, a small parameter e 
is related to a driving parameter: e — a 1 ! 3 . In this way the nonlinear equation 
under the zero initial data is derived for the leading order amplitude 

2iA'-2$'A + 37|A| 2 A = /, A(r)| r=0 = 0; (r = a 2/3 t). (5) 

A modulation of the complex amplitude A = A(t) in slow scale r = a 2 / 3 t is 
described by this equation. Our main discovery is a class of data /, under 
which a solution A(r) is infinitely increasing as r — > oo. This slow increasing is 
interpreted as the initial stage of the autoresonance. 



3 Rigid driving mode 

The phase function (3) at A = is considered in the section. Below this case 
will be referred to as rigid driving mode. 

It is expedient to represent a complex amplitude A by means of a pare of 
real functions R,ij)\ A = R(t) exp(i[^(r) — $(r)]). One can obtain explicit 
solution of the equation (5) for some special data /, 3>. For example, 

R(T) = \[f{C)dC, iP = t>-n/2, $( T ) = ?l£ R 2 (()dC (6) 

So in the case |/(r)| > fo > the amplitude increases linearly |i?(r)| = 
0{t), t — > oo if the driving phase $(t) is related with f(r) as implied by 
(6). As a result the leading order term of the asymptotic solution (4) has the 
order of unity u(t, a) — 0(1) at time r = 0(a~ 1//3 ), or what is the same, 
* = Oia- 1 ). 

In general case the existence theorem takes place 
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Theorem 1 Let the functions /(t),$(t) be continuous and the f(r) be uni- 
formly bounded for all r > 0. Then there exist a unique solution of the problem 
(6)forVr>0. 

This result gives a basis for construction of an asymptotics of the amplitude 
module R = R\t + i?o + 0(t -1 ' 2 ), t — > oo. Under assumption that the data 
have a smooth asymptotics at infinity: 

/(r) = /o+/it _1 +0(t- 2 ), $(t) = $ 3 T 3 +$ 2 T 2 +$ir+$o+0(r" 1 ), r oo. 

leading coefficients of an amplitude asymptotics are calculated: jR 2 = 2$3, 
3^fR\Ro = 2$2- So the autoresonance phenomenon, when Ri ^ 0, occurs under 
the conditions: 

7$ 3 >0, Jo 2 >8|* 3 |- (8) 

Here the second inequality is just the threshold condition for entering the au- 
toresonance. 

The first correction u\ in (4) is calculated from linearized equations. There is 
a secular term in m, which increases at infinity u\ — 0(t 3 ), t — > oo. Hence the 
first correction in expansion (4) remains lesser than the leading term a|tti| -C 
a 1 / 3 ^! just only over a time interval, which is not too long < t -C oT 1 . 

The main result in the rigid driving mode case is as follows 

Theorem 2 Let the right hand side in the equation (1) as given by (3) under 
A = 0, e = a 2 / 3 , (t = a 2 / 3 t) have the properties (7), (8). Then the system 
enters autoresonance. The amplitude of the leading order term in (4) found from 
equation (5) increases linearly |A(t)| = y / 2|$3|r + 0(1), r — > oo. Asymptotics 
(4) is available over a time interval, which is not too long < t < a~ 1+l/ , W > 
0. 

4 Loose driving mode 

The amplitude equation (5) is nonautonomous and in general case it can not be 
solved in the explicit form. However if both the coefficient $'(<5r) and the right 
hand side f(Sr) depend on the slow time St with a small parameter < 8 -C 1 
one can use an asymptotic approach known as adiabatic approximation. This 
case will be referred to as loose driving mode below and will be considered in 
this section. 

We consider the problem (1),(2) with phase driving function <p given by 
(3) with a small parameter e — a 2 / 3 and exponent A > 0; so the adiabatic 
parameter in amplitude equation (5) may be taken 8 = a 2X l 3 . Our main result 
is an asymptotics of the amplitude A = A(t, 8) as 8 — > 0, which is valid for 
large time r 3> S^ 1 . Using this asymptotics we find the condition under which 
the amplitude modulo |A(r, <5)| is infinitely increasing as St — > oo. This slow 
increasing is interpreted as the initial stage of the autoresonance. 
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Let us go over to the asymptotic constructions. First, the scaling transfor- 
mation is performed in the equation (5) 

A(r,S) = (f/3 1 )^B(( 1 S), C = V/S, V = fj d£. (9) 

Then the problem for the complex amplitude B = P + iQ is represented by a 
system of differential equations for two functions, namely for real and imaginary 
parts 

Q' -[P 2 + Q 2 -n( V )}P = -l + SF( v )Q, Q(£*)|c=o = 0, (10) 
P' + [P 2 + Q 2 - Q( V )}Q = SF( V )P, P(C,%=o = 0. (11) 

Coefficients 

fifa) = 2$'/(3 7 ) 1/3 / 2/3 (^), F( V ) = -2/7(3 7 ) 1/3 / 5/3 (5r) 

depend on the new slow time r\ = < 5 C 1. These equations are considered 
either on the semiaxis ( > or ( < depending upon the signs of nonlinearity 
7 = ±1. 

A small parameter S occurs as the factor in equations (10), (11) and it defines 
the slow time in data as well. To solve such perturbed problem we apply an 
asymptotic method, which is ascribed to Bogolubov, Krylov, Kuzmak, Haber- 
man, [8,11,12]. The basis of the approach is a two-parametric periodic solution 
of the unperturbed problem under frozen data fl = const. In the case under 
consideration the unperturbed system (as 6 = 0, O = const) is autonomous and 
hamiltonian: 

Q' - [P 2 + Q 2 - Sl]P = -1, 

P' + [P 2 + Q 2 - Q]Q = 0. 

The system is integrable and an explicit integral representation of the solution 
can be obtained. But there is no need in an explicit formula to grasp the 
situation and to find the condition for entering the autoresonance. To this end 
the Hamiltonian as the first integral may be used 

j(P 2 +Q 2 -nf - P = const. 

One can see the family of periodic solutions from the picture of the phase space 
portrait, where nearly all trajectories are closed except for two embedded sepa- 
ratrixs. To parametrize the solution it is reasonable to use both the plane area II 
covered by trajectory and the phase shift s in the form P,Q = Po, Qo(C+s, II, SI). 
Moreover the solution depends on the parameter Jl since O is present in the 
equations. 

The behavior of separatrixs depending upon the Jl has crucial role in detec- 
tion of the autoresonance. Both these lines tend to the large circle P 2 + Q 2 = 
as Q — ► oo, while the distance between them is small and has the order 
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C(ri -1 / 4 ) outside some neighbourhood of the saddle point. However, the area 
between separatrix loops increases and it is minorized by a large magnitude 
mfl 1 / 4 , (to = const > 0), since the separatrix length has the order of ©(fi 1 / 2 ). 
In the considered case with frozen data the solution Pq , Qq starting from zero 
is periodic and finite at any time (. Its phase space trajectory depends on the 
parameter f2. It may be located either inside the inner separatrix or between 
them. 

We return to the perturbed problem (10), (11) in which either Q(r]) ^ const 
and 6F(rj) ^ or only one of the inequalities f2(ry) ^ const, 8F(rf) ^ holds. 
One can construct an asymptotic solution which is valid over a long time interval 
by using both the adiabaticity of the coefficient SI and the smallncss of the factor 
S. Following [7-9] we introduce the two-scale representation of the solution. In 
this approach the leading order term is taken as the unperturbed solution with 
additional dependence on the slow time 

Qo = Qo{<r, 11(77), SI( V )), P = Po(v, nfa), 0(r?)) 

while the fast variable is a = S Si(rj) + So (77). So an asymptotic solution of 
(10), (11) is represented by a pare of functions which are periodic with respect to 
fast variable a and depend on slow time by means of II, Sl(rj). In this approach 
it is necessary to find slow varying functions II, S\, So(r]) in order to identify the 
leading order term of the asymptotic solution . 

We are interested in the solution in which the amplitude \P\ + \Q\ is an 
increasing function of the slow time 77. Evidently, the structure of the fast 
variable a is not needed at this stage. The main result in the section is based on 
finding the slow varying function 11(77). K is obtained from differential equation, 
[11,12] which turns up trivial due to the appropriate parametrization 

a*n = F(0)n, n(7?)| tf=0 = n . 

Here n is the area covered by the initial trajectory (as r\ = 0). In this way a 
remarkably simple formula is obtained: 

H(77)=H |/(0)//(5t)[. 

In the special case when / = const the area is an adiabatic invariant n = const. 

Thus using (6), one can write down an adiabatic approximation of the solu- 
tion of the amplitude problem (5) as follows 

A(t,S) = (f(5r)/3 1 ) 1 / 3 [Po( f yMv),^(v))+^Qo('jMv)M(v))} + Mr,S), (12) 

Here the remainder A(r, S) — > as S — > 0. As to the remainder it is known 
[13] that it is evaluated through small parameter A(r, 5) = 0{8), 8 — > and 
the estimate is uniform over long time interval < ( < 770 {Vo = const > 
0). But this strict result is not sufficient for us, because the autoresonance 
phenomenon is detected only at very far time ( 3> S^ 1 where the leading term of 
the asymptotics (12) may be increasing. There is no conflict since the remainder 
is small A(r, 5) = 0(1), S — > over a longer interval, on which the slow time 
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can become large 77 = 5£ ^> 1. However, one has to keep in mind that the 
remainder becomes worse with the growth of time. Hence the approximation 
given by (12) is valid so far as the remainder can be neglected with respect to 
leading term. Analysis of the first correction in the adiabatic asymptotic may 
suggest a limiting time interval but we do not treat corrections in this work. 
One can say a priori that an adiabatic approximation of the amplitude is valid 
until the small amplitude expansion (4) is suitable. 

The main result of the section is derived from the formula (12) under the 
relation S = a», (fj, = 2 A/3). 

Theorem 3 Let the right hand side in the equation (1) have both the amplitude 
and the phase f = f(6), ip = t + a~ (9 = a» +2 / 3 t, (j, > 0) which satisfy 
1/(0)1 > |/(0)|; I*' (0)| -> 00 as 9 00; let the function /7 1 / 3 / 2 / 3 (0) -» 00 
be monotonous and increase to infinity of 9. Then entering the autoresonance 
depends on the initial value of the parameter f2(0) = 2<E>'(0)/(37) 1 / 3 / 2 / 3 (0). 
Autoresonance does not occurs if 17(0) > 3/2 1 ' 3 . Autoresonance occurs under 
fi(0) < 3/2 1 / 3 . In the last case the leading order term of an amplitude asymp- 
totics is determined by the driving frequency as follows 

|A(t,o")| = (^'(e)|) 1/2 + 0(|$T 1/4 ) + OK), <* -0, 9^oo. 

In order to prove the theorem it is enough to observe a slow deformation of 
the phase plane trajectory Pq, Qo{c, U(rf), Cl(rf)) under variation of the parame- 
ter 77. In the first case the trajectory is located inside loop of the inner separatrix 
and it encircles the equilibrium point near zero. In the second case the trajectory 
is located between separatrix loops and it encircles another equilibrium point 
near P = V^M- Under given conditions the plane area n covered by trajec- 
tory is not increasing, while separatrix loops grows like the circle P 2 + Q 2 = Q 
as ft — > 00. Hence in the first case the trajectory remains in bounded part of 
the phase plane for all time. In the second case one has to take into account 
that the plane area between loops is increasing. Hence the trajectory remains 
between loops at any time and it is going to infinity as fl — > 00. 
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